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Abstract. We examine the results linking cosmic ray flux (CRF) variations to global 
climate change. We then proceed to study various periods over which there are estimates 
for the radiative forcing, temperature change and CRF variations relative to today. These 
include the Phanerozoic as a whole, the Cretaceous, the Eocene, the Last Glacial Max- 
imum, the 20* century, as well as the 11-yr solar cycle. This enables us to place quan- 
titative limits on climate sensitivity to both changes in the CRF, ^cRi a-nd the radia- 
tive budget, F, under equilibrium. Under the assumption that the CRF is indeed a cli- 
mate driver, we find that the sensitivity to CRF variations is consistently fitted with fi = 
-^o{dTgiobai / d^CB.) ~ 6.5±2.5°K (where <i>o is the CR energy flux today). Addition- 
ally, the sensitivity to radiative forcing changes is A = dTgiobai/dF\^^ — 0.35±0.09°KW^^ m^, 
at the current temperature, while its temperature derivative is negligible with {dX/dT)o = 
0.01±0.03m^ W"^. If the observed CRF/climate link is ignored, the best sensitivity ob- 
tained is A = 0.54±0.12°KW-im2 and (dA/dT)o = -0.02±0.05m2 W'^. The CRF/climate 
link therefore implies that the increased solar luminosity and reduced CRF over the pre- 
vious century should have contributed a warming of 0.37±0.13°K, while the rest should 
be mainly attributed to anthropogenic causes. Without any effect of cosmic rays, the in- 
crease in solar luminosity would correspond to an increased temperature of 0.16±0.04°K. 



1. Introduction 

Accumulating evidence suggests that solar activity is re- 
sponsible for at least some climatic variability. These in- 
clude correlations between solar activity and either direct 
climatic variables or indirect climate proxies over time scales 
ranging from days to millennia [Herschel, 1796; Eddy, 1976; 
Labitzke and van Loon, 1992; Friis-Christensen and Lassen, 
1991; Soon et al, 1996a, 2000; Beer et al, 2000; Hodell et at, 
2001; NeJJ et al, 2001]. It is therefore difficult at this point 
to argue against the existence of any causal link between 
solar activity and climate on Earth. However, the climatic 
variability attributable to solar activity is larger than could 
be expected from the typical 0.1% changes in the solar ir- 
radiance observed over the decadal to centennial time scale 
[Beer et al, 2000; Soon et al., 2000]. Thus, an amplifier is 
required unless the sensitivity to changes in the radiative 
forcing is uncomfortably high. 

The first suggestion for an amplifier of solar activity was 
suggested by Ney [1959], who pointed out that if climate is 
sensitive to the amount of tropospheric ionization, it would 
also be sensitive to solar activity since the solar wind mod- 
ulates the cosmic ray fiux (CRF), and with it, the amount 
of tropospheric ionization [Ney, 1959]. 

Over the solar cycle, the solar wind strength varies con- 
siderably, such that the amount of tropospheric ionization 
changes by typically 5%-10%. Svensmark [1998, 2000], 
Marsh and Svensmark [2000a] as well as Palle Bago and 
Butler [2000] have shown that the variations in the amount 
of low altitude cloud cover (LACC) nicely correlate with 
the CRF reaching Earth over two decades. Both signals, of 
cloud cover and CRF variations lag by typically half a year 
behind other solar activity indices. This suggests that it is 
more likely that the cloud cover is directly related to the 
CRF than directly to solar activity. 
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More recent data on the LACC seems to exhibit a weaker 
correlation with the variable CRF (e.g., Farrar [2000]). 
There are however a few peculiarities in the data which are 
indicative of a calibration problem, which once removed, 
seem to recover the high correlation between the CRF and 
the LACC [Marsh and Svensmark, 2003]. For an objec- 
tive review, the reader is encouraged to read Carslaw et al. 
[2002]. 

The above correlations between CRF variability and cli- 
mate (and in particular, cloud cover), indicate that CRF 
modulations appear to be responsible for climate variabil- 
ity, most probably through modulation of the amount of 
LACC. Nevertheless, since all of the above CRF variability 
ultimately originates from solar activity changes, it is not 
possible to unequivocally rule out the possibility that the 
CRF/climate correlations are coincidental, and that both 
are independently modulated by solar activity with similar 
lags. 

An independent CRF/climate correlation on a much 
longer time scale, in which variations in the CRF do not orig- 
inate from solar variability, was found by Shaviv [2002a, b]; 
Shaviv and Veizer [2003]. It was shown using astronomical 
data that a large ~ 0(1) CRF variability should arise from 
our passages through the galactic spiral arms, with a period 
of 132 ± 25 Ma. It was also shown that the CRF history 
can actually be reconstructed using the cosmic-ray expo- 
sure age data of Iron meteorites, exhibiting a periodicity of 
143 ± 10 Ma and a phase consistent with the astronomical 
data. Moreover, it was found that the reconstructed CRF 
nicely synchronizes to the occurrence of ice-age epochs on 
Earth, which appeared on average every 145 ± 7 Ma over 
the past billion years. Additionally, the mid-point of the 
ice-age epochs is predicted to lag by 31 ± 8 Ma after the 
mid-point of the spiral arm crossing, while it is observed to 
lag by 33 ± 20 Ma. That is, the CRF and ice-age epoch 
signals agree in both phase and period. The same analysis 
also revealed that the long term star-formation activity of 
the Milky Way correlates with long term glacial activity on 
Earth. In particular, a dearth in star formation between 1 
and 2 Ga before present, coincides with a long period during 
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which glaciations appear to have been totally absent [Sha- 
viv, 2002b, 2003]. 

We should also point out several experimental results sup- 
porting, though not proving yet, a CRF/cloud cover link. 
Harrison and ApUn [2001] found experimentally that CN 
formation is correlated with natural Poisson variability in 
cosmic ray showers. In other words, this link appears to be 
more than hypothetical. In another set of experiments, it 
was shown that cosmic rays play a decisive role in the forma- 
tion of small clusters [Eichkorn et al, 2003]. If these small 
clusters can bo shown to grow quickly enough, as opposed 
to being scavenged by large particles, the link between cos- 
mic rays and the formation of cloud condensation rmclei and 
ultimately cloud cover could be firmly established. 

We will not dwell here on the actual mechanism responsi- 
ble for CRF link with cloud behavior. We will simply assume 
henceforth that this link exists, as supported by empirical 
and experimental data, even though it is still an issue of 
debate. This point has to be kept in mind since the conclu- 
sions we shall reach, will only be valid if this assumption is 
correct. 

Using the above assumption, we study several time 
scales to see whether estimates on global temperature sen- 
sitivity can be placed, together with estimates on the 
CRF/temperature relation. We will do so by comparing 
the observed temperature changes with changes in the ra- 
diative budget, an approach previously pursued in numerous 
analyses (e.g., Hoffert and Covey [1992]; Covey et al. [1996]; 
Hansen et al. [1993]; Cregory et al. [2002]). This method 
for obtaining the global temperature sensitivity using palco- 
data is orthogonal to the usage of global circulation models 
(GCMs) upon which often quoted results are based [IPCC, 
2001]. Hence, the two methods suffer from altogether differ- 
ent errors. It is therefore clearly advantageous to follow this 
path as an independent estimate. For example. Cess et al. 
[1989] have shown that the large uncertainty in the sensi- 
tivity obtained in GCMs stems from the uncertain feedback 
of cloud cover. Since we use the actual global data, all the 
feedbacks are implicitly considered. The main contribution 
in this work is to specifically consider the contribution of the 
CRF to the changed radiative budget. As a note of caution, 
one should keep in mind that the most notable assumption 
in this method is the quantification of climate sensitivity 
with one number. In other words, it assumes that on av- 
erage Earth's climate responds the same irrespective of the 
geographic, temporal or frequency space distribution of the 
radiation budget changes. It also assumes that different ra- 
diative forcings act linearly. 

Once the radiative forcing and temperature changes are 
obtained, the sensitivities can be estimated with 
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AF, which is the globally averaged change in the radiation 
fiux (per unit surface axea), will also include here the con- 
tribution AFcRF arising from a changed energy flux $cfl 
of cosmic rays. Note also that over short time scales, AT or 
AF have to be properly modifled to include the finite heat 
capacity of the system, and the consequent finite adjustment 
time it has. We should also consider the possibility that A 
is dependent on the temperature. For example, the positive 
climate feedback arising from the formation of ice sheets 
could increase the sensitivity of a glaciated Earth, while the 
reduced atmospheric water vapor content, can reduce the 
sensitivity. 

In addition to A, we will also estimate the sensitivity to 
CRF variations, defined as: 
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where is the cosmic ray energy flux reaching Earth at 
energies of ~ 10 GeV (the energies responsible for tropo- 
spheric ionization), while $o is the average flux reaching 
today. 



2. Radiative forcing of low altitude cloud 
cover 

Without a detailed physical model for the effects of cos- 
mic rays on clouds or a detailed enough record of radiation 
budget meeisurements correlated with the solar cycle, it is 
hard to accurately determine the quantitative link between 
CRF variations and changes in the global radiation budget. 
In particular, it is hard to do so without limiting ourselves 
to various approximations. Nevertheless, this link is impor- 
tant since it will be used in most of our estimates for the 
global temperature sensitivity. 

The basic observation we use to estimate the radiative 
forcing of clouds is the apparent correlation between CRF 
variations and the amount of low altitude cloud cover. A 
naive approximation is to assume that the whole climatic ef- 
fect can be described by variations in the extent of the cloud 
cover, namely, that we neglect effects in the cloud proper- 
ties, or possible climatic effects associated with atmospheric 
ionization but not with clouds. It also implies that the geo- 
graphical distribution of the effect is the same as low altitude 
clouds on average. We will first estimate this "zeroth" order 
term and then try to estimate the possible contribution of 
other corrections. 

Amount of cloud cover: Over the solar cycle, the varying 
CRF appears to cause a 1.2 to 2.0% (absolute) change in the 
amount of LACC [Marsh and Svensmark, 2000b; Ktrkby and 
Laaksonen, 2000; Marsden and Lingenfelter, 2003; Carslaw 
et al., 2002]. We will therefore adopt a change of 1.6 ±0.4% 
in the LACC. 

The total radiative forcing of the LACC is estimated to 
be — 16.7Wm~^ from the average 26.6% cloud cover [?]. If 
one however compares the forcing of the total cloud cover 
from different hemispheres and different experiments (Nim- 
bus and ERBE, Ardanuy et al. [1991]), one finds variations 
which are typically 2.5Wm^^ on the ~ SOWm^'^ short- 
wave (SW) "cooling" and 7Wm^'^ on the 27Wm"2 long- 
wave (LW) "warming" . Since LACC typically comprise half 
of the total amount cloud cover, an error of ~ 4Wm~^ is to 
be expected. 

Thus, the changed radiative forcing AFf associated with 
the varying amount of cloud cover, should be — 1.0 ± 
0.35Wm"^. This implicitly assumes that the incremental 
cloud cover has the same average net radiative properties as 
the whole 27% of the LACC. 

Cloud Optical Depth: Changes in the cloud properties 
could take place in addition to changes in the cloud amount. 
According to Marsden and Lingenfelter [2003], there is a 
small negative correlation between the average LACC opac- 
ity r and the varying CRF. Over the solar cycle, f changes 
by —4% relative to its global average of about 4 in regions 
defined to be covered by LACC. 

Is such a change in f reasonable? According Marsden 
and Lingenfelter [2003], there are two limiting cases for 
the effects on cloud properties. The first is changing the 
number density of cloud condensation nuclei (CCN) given 
a fixed amount of Liquid Water Content (LWC), that is, 
CCN limited. This is similar to the "Twomey effect" where 
enhanced aerosol density affects the droplet size and cloud 
albedo [Twomey, 1977; Rosenfeld, 2000]. The second case 
is increasing the CCN density together with the LWC, and 
obtaining similar sized drops (LWC limited) . Although the 
two cases are plausible, they do not change the cloud prop- 
erties in the same way. 
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One can show that a cloud's optical depth for SW ab- 
sorption is (e.g., Marsden and Lingenfelter [2003]) 



3 PobAz 

2 poRcB ' 



(3) 



where po is the density of water, p^s is the mass load- 
ing of water (i.e., its effective density), Az is the vertical 
extent of the cloud and R^g is the efitective radius of the 
cloud droplets, defined as the ratio between the 3'^'' and 2"'' 
moments of the droplet distribution {{r'^') / (^'^))- In the 

case of a CCN limited condensation, Rcs cx ^(^QjY , and t 
will increase with nccN, while in the LWC limiting case, 
Peff oc nccN and t will increase as well. One can therefore 
write: 



5r _ ^ SnccN 
T nccN 



(4) 



with /3 = 1 for LWC limited case and = 1 /3 for the CCN 
limited case. Since the lower troposphere ionization rate 
changes by about 7% between solar minimum and maxi- 
mum, we should expect to get at most a similar increase in 
the CCN. Thus we should expect (5r/r^7%. The fact that a 
negative change in r was observed [Marsden and Lingenfel- 
ter, 2003], could arise because the increase in cloud lifetime 
results with thinner clouds on average. 

Next, one can approximate the relation between t and 
cloud albedo A, by the relation [Hobbs, 1993]: 
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where ri/2 « 6.7 for an asymmetry parameter of 0.85 
[Hobbs, 1993]. Once we differentiate, we find: 
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If wc consider the transmission T ~ 0.75 of the atmosphere 
(to obtain a top-of-atmosphcrc albedo, from a cloud-top 
albedo), that the LACC covers only a fraction fiow of the 
globe, and that the average top-of-atmosphere incidence of 
radiation is F = 344Wm^^, wo find that the change in 
albedo is responsible for a changed radiation budget of: 



^-^FT'^fioyjSf « -K0.13Wm-^ 



(7) 



The positive sign implies that the small apparent reduction 
in f contributes a small warming contribution. 

If we had no knowledge of r, changes in it could have 
resulted with a correction to AF^ which are only as large 
as -0.23Wm-2 (for the LWC limited case, and 5r/r<7%). 
We take this uncertainly in f as another source of error for 
the radiative forcing AF. 

Cloud Emissivity: There could still be more physical 
terms contnbutmg to AF. If the LWC in the clouds can vary 
as well (that is, the clouds are not CCN limited but rather 
water limited), then also the IR emissivity can change. It 
will do so by changing the emissivity, relative to black body 
(e.g., Stephens [1978]) which is given by: 



€{6z) = 1 - exp(-r/ij). 



(8) 



where we have defined tir = aoPesSz. Here, peff is the liq- 
uid water content, ao is the mass absorption coefficient (for 
water clouds, ao w 0.13m^g~^ [Stephens, 1978]) and Sz is 
the width of the cloud layer above a given point. By chang- 
ing the emissivity, we change the outgoing long-wavelength 
flux by 

AFiH = TaT^fio^Ae = eM-riR)TaT\iR^^ (9) 

Pefl 



where T ^ 0.6 is the transmittance of the atmo- 
sphere to IR, above the cloud. For typically small val- 
ues of Peff of 0.3g/m^, and 5z = 100m (which would 
give the largest effect), we get corrections of AFir = 
l.lWm"^(Apeff/pefE)^0.1Wm"^. This positive flux out- 
wards tends to cool (i.e., increase the CRF/temperature ef- 
fect), but it is a small effect. 

By changing the emissivity, wc can also shift the appar- 
ent location of the top of the clouds, and with it their tem- 
perature. In other words, we should expect outgoing LW 
radiation to come from higher up the atmosphere where the 
temperature is lower. 

A higher petF will shift the IR emission "surface" vertically 
by typically: 
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Using the black body law, the change in the radiative emis- 
sion over the solar cycle will therefore be less than 

-AFAT^/w4aT3A^^^0.2^Wm-^ (11) 

az Peff 

once globally averaged. For the last inequality, wc took a 
typically low Peff of 0.3g/m^, a wet adiabat of dT/dz ~ 
0.6°is:(100m)"^ and fio-w ~ 0.28. Since Apefi/peff <0.1, this 
effect will be even smaller at best (and in opposite sign as 
the previous effect). 

This result is also reasonable considering that the total 
long wavelength heating effect of LACC was estimated to 
be ~ 3.5Wm^'^ [?], while cloud albedo is responsible for 
a globally averaged cooling of ^ 20Wm^^, implying that 
changes in albedo will likely be more important for chang- 
ing the radiative budget arising from LACC vaxiations. 

Ocean Bias: Additional unaccounted effects are possible. 
For example, as we have previously stated, we implicitly as- 
sume that the incremental change in the LACC contributes 
to the radiative budget as the LACC on average. This need 
not be the case. Suppose the CRF-LACC effect primarily 
takes place over the oceans. This is reasonable because it 
is the oceans where the density of CCN is lowest. Since the 
oceans have a lower albedo than land, covering ocean sur- 
face has a larger net decrease in the radiative budget than 
covering land mass. 

If the land albedo is typically higher by 15% and LACC is 
typically half as frequent over land, then the albedo change 
arising from uncovering clouds over only oceans is higher 
than the albedo change obtained when uncovering the av- 
erage LACC, by A^ ^ 3%. This corresponds to a flux 
change of AF^a ~ AAT'^AfiowF « -hO.lWm'^ That is, 
a likely ocean bias implies that we are slightly underestimat- 
ing AFcflF. 

Other effects: If the effect is geographically localized to 
certain areas, then a larger discrepancy could arise if the 
radiative properties of the LACC over those geographic re- 
gions is signiflcantly different from the properties of LACC 
on average. A correlation map between LACC variations 
and CRF change [Marsh and Svensmark, 2000b], reveals 
that some regions (particularly over oceans) stand out with 
a higher correlation than others. Nevertheless, they do not 
appear to cluster around particular latitudes or other spe- 
cial regions. Thus, the assumption of geographic uniformity 
may bo not that bad. 

Another hard to estimate effect could arise from the ex- 
pected increase in cloud lifetime. For example, cumulus- 
type clouds could penetrate into higher altitudes, thereby 
reducing their IR emission. 



X- 4 



SHAVIV: COSMIC RAYS AND CLIMATE SENSITIVITY 



Thus, until wc fully understand all the details in the 
physical picture, we should take the estimated radiative 
forcing and the error with a grain of salt. Taking the 
above into considerations, our best estimate for the radia- 
tive forcing of the cloud cover variations over the solar cy- 
cle is AFcRF = AFf + AFa + AFaa = 1-0 ± 0.4Wm~^ 
globally averaged (we have neglected AF/ji<0.1 Wm~^ and 
-AFAT<0.02Wm~^). This should be compared with the 
0.1% solar flux variations, giving rise to an extra "direct" 
forcing of AFaux w 0.35Wm"^ [Prohlich and Lean, 1998]. 

A related number describes the relation between CRF 
variations and changes in the low altitude ionization rate. 
At energies which can reach the lower troposphere and low 
geomagnetic latitudes (>10G'eV^), a fair representative of 
the flux are the neutron measurements at Hunacayo Peru 
and Halaekala, Hawaii (near the geomagnetic equator, with 
a rigidity cutoff of 12.9 GcV). There, the amplitude in the 
CRF variations over the past few solar cycles was 7 ± 2% 
[e.g., Svensmark, 1998] Thus, we find that the radiative sen- 
sitivity to CRF variations is about 



dF 

' d^CR 



14±6Wm" 



(12) 



3. Estimating Climate Sensitivity 

We now proceed to estimate the climate sensitivity. We 
do so by comparing the radiative forcing change between two 
eras to the temperature change which ensued, using eq. 1. 
In most estimates, we will rely on the results of §2 to ob- 
tained the contribution of the changed CRF to the changed 
radiative forcing. These include seven different comparisons, 
spanning from variations over the solar cycles, to variations 
over the Phanerozoic as a whole. Subsequently, we will com- 
bine the results to obtain our best estimate for the climate 
sensitivity. 

3.1. The AT/CO2 correlation over the Phanerozoic 

Shaviv and Veizer [2003] have shown that more than 2/3's 
of the variance in the reconstructed tropical temperature 
variability ATi^op over the Phanerozoic can be explained 
using the variable CRF, which could be reconstructed using 
Iron meteorites. On the other hand, it was shown that the 
reconstructed atmospheric CO2 variations do not appear to 
have any clear correlation with the reconstructed tempera- 
ture. The large correlation between reconstructed CRF and 
temperature is seen in fig. 1. It is this correlation which led 
the authors to conclude that the Phanerozoic climate is pri- 
marily driven by a celestial driver. The lack of any apparent 
correlation with CO2 was used to place a limit on the global 
climate sensitivity. 

A subsequent analysis by Royer et al. [2004] has shown 
that pH corrections could have been important at offsetting 
the (5^*0 record upon which the temperature reconstruction 
is based. In particular. The pH correction term of Royer 
et al. [2004] has the form: 

ATj,H = a{log Rco2 + log K{t) - log n{t)}, (13) 

where Rco2 is the atmospheric partial pressure of relative 
to today, A(f) is (Ca)(i)/(Ca)(0) - the mean concentration 
of dissolved calcium in the water relative to today, while 
Q.{t) is [Ca'^'^][CO^~]/ Ksp at time t relative to today. The 
value of a obtained in Royer et al. [2004] assumes no ice- 
volume correction, yielding a = 3.4°K. Once the ice-volume 
effect on 5^^0 is considered [Veizer et al., 2000; Shaviv and 
Veizer], one obtains: a « 1.4°K. 

Since the pH correction depends on Rco2, so will the cor- 
rected temperature. A simple correlation between the cor- 
rected temperature and the reconstructed Rco2 will then be 
meaningless. Instead, the method to proceed is to defined a 



CO2 "uncorrected" temperature as: 

AT' = AT - a\ogRco2. 



(14) 



Any correlation that this signal will have with Rco2 will 
then be real, since this "temperature" depends only on 5^*0, 
and the small A and f2 terms. This uncorrected temperature 
can then be fitted with 



A + Bt + {C -a logio 2) log^ Rco2 
+Dg{^CRF{t)). 



(15) 



A and B allow for systematic secular trends in the data. D 
relates the cosmic ray energy flux $(f) to AT [Shaviv and 
Veizer, 2003]. The term (alogjg 2)log2Rco2 was added such 
that C will keep its original meaning in Shaviv and Veizer 
[2003], which is the tropical temperature increase associated 
with a doubled Rco2. 

The lack of a correlation between 5^*0 and Rco2 [Shaviv 
and Veizer, 2003] originates from the fact that (C— alog^Q 2) 
happens to be coincident ally close to 0. In other words, 
the pH correction to 5^^0 and AT, happens to be similar 
to the tropical temperature sensitivity to changes in Rco2 
(Without the pH correction, the preferred value for C in the 
absence of correlation is not alog]^Q2, but 0). Scientifically, 
this is somewhat unfortunate, because without this coinci- 
dence the <5^**0 signal would have had a clear correlation 
with the Rco2 signal, and the Rco2 fingerprint would have 
boon discernible in the Phanerozoic data. 

If wc repeat the analysis of Shaviv and Veizer (2003), and 
consider also the effects of alog A(t) — log n(t) introduced 
by Royer et al. [2004], and corrected for Rco2 as described 
above, we obtain: C = 0.69° K (or an upper limit of 1.12, 
1.42 and 1.73°K at 68%, 90% and 99% confidence levels, 
and a lower limit of 0.39, 0.10, -0.21°K, respectively). This 
gives A = 0.28 ± 0.15°K W"^ m^. 



r=-0.77 
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Figure 1. The high correlation between the reconstructed 
temperature and CRF over the Phanerozoic can be used 
to estimate global sensitivity. Here, AT' is the recon- 
structed temperature AT of Veizer et al. [2000] over the 
past 550 Ma, with a small linear temperature increase 
of 1.7°K(t/550Ma) subtracted (see Shaviv and Veizer 
[2003]). The CRF is one of 3 reconstructions used in Shaviv 
and Veizer [2003]. The two others differ in the total am- 
plitude of variations. The two independent signals have a 
high Pearson correlation of r = —0.77. Although a statisti- 
cally significants limits on p cannot be placed, lower p's are 
favored (with p ~ 0.3 producing the best fit). Nevertheless, 
the value of p = 0.5 is theoretically preferred. The data is 
smoothed with a 50 Ma moving average. Points are 10 Ma 
intervals. 
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Without the effect of cosmic rays (i.e., witfi tlic D term 
removed in the model given by cq. 15), more of the re- 
constructed temperature variabihty can be explained with 
CO2, and the estimate range for A broadens respectively to 
A = 0.36±0.22°KW"^m^. More limits are given in table 1. 

Note that this estimate is independent of a determined in 
§2. The first range for A simply assumes that a CRF/climate 
link exists, while the second quoted range, even neglects this 
assumption. 

3.2. CRF/AT correlation over the Phanerozoic 

The significant correlation between CRF and tempera- 
ture over the Phanerozoic was also used to place limits on 
the ratio between CRF variations and temperature change. 
Together with the results of §2 we can place a limit on A. 

In Shaviv and Veizer [2003], it was found that if ATtrop 
is approximated with ATirop = D [{'^cr/^oY ~ 1]) where 
p = 1/2 and $0 is the CRF today, then D = 8± 4°K. Al- 
most all the contribution to the error arises from our limited 
knowledge of the actual amplitude of the CRF variations. If 
we generalize to other power laws p between 0.25 to 1.5, and 
repeat the procedure described in Shaviv and Veizer [2003] , 
we find that 



dTglobal 



d^c 



= pD-d, 

*CR=*0 01 trap 



= 6.0±3.0°K, (16) 



where we have taken STgiotai / STtrop ~ 1-5 as typically ob- 
tained in GCMs [IPCC, 2001]. The reason the error docs 
not increase much once wc introduce a range of p's is be- 
cause fi (but not D) is rather insensitive to p. Note also 
that p = 1/2 is theoretically preferred [Yu, 2002; Harrison 
and Aplin, 2001; Ermakov et al, 1997], and a small p is also 
favored by the empirical data (sec fig. 1). 

Using the result for a obtained in §2, we find A = n/a = 
0.441q;2| KW~^m^. where we quote the median A and the 
le**" and 84*'' percentiles (l-a). More details on the distri- 
bution appear in fig. 3 and table 1. 

3.3. Bounds from the total T and ^CR variations 
over the Phanerozoic 

Using the same Phanerozoic data and an altogether dif- 
ferent set of argumentations, we can place additional limits 
on A and on p,. We do not know accurately how large are 
the absolute CRF variations that give rise to the temper- 
ature oscillation over the Phanerozoic (hence the relatively 
large error on fj, and D). Nevertheless, we know that there 
is a maximum increase of ~ 2°K in the tropical temperature 
above today's tropical temperature, once averaged over the 
50 Ma time scale [Veizer et al., 2000]. This approximately 
corresponds to an increase of AT ~ 3°K globally. We as- 
sume that this could arise by removing at most fmin ~ 80% 
of the LACC, which would give rise to a global cooling of 
fmin^FLACc « /minlTWm"^. Namely, 



An 



AT 



' fminAFhACC 



0.22°KW 



-1 2 
m . 



(17) 



This is an absolute minimum for the climate sensitivity. 
Otherwise, the CRF-temperature link observed over the 
Phanerozoic would require too large a radiative budget 
change to be explained by LACC variations. 

A similar limit can be placed on fx. The largest in- 
crease over the Phanerozoic of AT ~ 3°K relative to to- 
day, and which can be attributed to CRF variation, can- 
not arise from a flux decrease larger than 85%, since larger 
reductions are inconsistent with various astrophysical con- 
straints. This implies /i>3.0p/(l — 0.15^)°K, if a relation of 
the form AT oc exists. Likewise, the maximum temper- 
ature decrease over the Phanerozoic is about 2.5°K over the 
tropics, or ~ 3.75° K globally, while the astronomical con- 
straints give a maximum flux increase of about 95% relative 



to today for models still satisfying astronomical constraints 
[Shaviv, 2002b]. This gives /i^3.75p/(1.95^ - 1)°K. Com- 
bining the two, wc get /i>3.7°K if p is unconstrained, or 
^i>4.7( or 4.0)°K if p = l/2( or 1). 

Using §3.2-3.3, we thus obtain /i = 6.5 ± 2.5°K. 

3.4. Cretaceous and Eocene Climates 

Particular geological epochs were studied under more 
scrutiny, and without being averaged out on the 50 Ma time 
scale, as the Phanerozoic data was. In particular, there 
were estimates for both the radiative forcing and temper- 
ature change of two geological periods which were partic- 
ularly warm relative to today. One is the Cretaceous, at 
about 100 Ma before present, and the second is the Eocene, 
at 50 Ma. 

Barron et al, [1995] estimated that the mid-Cretaceous 
was 7 ± 2°K warmer than today, and that it arose from an 
increase of 8 ± 3.5 Wm~^ in the radiative budget. However, 
Covey et al. [1996] point out that this estimate included 
only the change from the increased amount of atmospheric 
CO2 and it did not include the increased forcing associated 
with surface albedo changes. Once taken into account, the 
Hoffert and Covey [1992] estimate for the radiative forcing 
increases to 15.7 ± 6.8Wm~^. Their estimate for the tem- 
perature increase is also larger at 9 ± 2°K. We adopt the 
Hoffert and Covey [1992] estimate here. 

Covey et al. [1996] estimated the temperature and radia- 
tive flux increases associated with the Eocene. They are 
4.3 ± 2.1°K and 11.8 ± 3.6Wm''^ respectively. The tem- 
perature and forcings can be used to estimate the sensi- 
tivity through A — AT/AT. The results for the Eocene 



and Cretaceous are A = 0.37' 



+0.20° 



KW" 



and A = 



0.57l2;t8°KW-im2 respectively. They are also summarized 
in fig. 2 and table 1. 
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Figure 2. The Probability Distribution function for A (and 
ATx2) obtained by comparing radiation budget differences 
to temperature change over various time scales, assuming 
that CRF variations do not affect the global climate (though 
it does includes the small solar luminosity changes). We 
also assume, as Gregory et al. [2002], that AT and AF 
entering A have Gaussian errors. The cases are: (1) Tem- 
perature increase over the past century (following Gregory 
et al. [2002]). (2) Temperature variations over 300 yrs of 
solar cycles. (3) Warming since the LGM (following Hoffert 
and Covey [1992] and Hansen et al. [1993]). (4) Cooling 
relative to the Cretaceous (~ 100 Ma, Hoffert and Covey 
[1992]). (5) Cooling relative to the Eocene (~ 55 Ma, 
following Covey et al. [1996]). And (6) Phanerozoic AT 
vs. Rco2 (§3.1, following Shaviv and Veizer [2003]). We 
assume that the temperature increase ATx2 following the 
doubling of the atmospheric CO2 content corresponds to an 
increase of 3.71 Wm'^ [Myhre et al., 1998]. 
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These estimates do not include however the possibihty 
that CRF variations affect chniate. To estimate this effect, 
we estimate the CRF differences between the two geological 
periods and today using the CRF reconstruction described 
in Shaviv [2002b] . We then calculate AFcr arising from the 
CRF change and use AT = A(AFo + AFcr) to obtain A. 

We find that the CRF was 20% to 70% of the flux today 
during the mid-Cretaceous. Through the effects on clouds, 
this should have contributed towards an increase in temper- 
ature, and therefore reduce our estimate for A. During the 
Eocene, $ch should have been between 0% and 20% higher 
than today. From the 6 epochs described here, this is the 
ordy case in which the effect of the CRF is to increase the 
estimate for the cliinatc sensitivity. 

The radiative forcing associated with the CRF variations 
can be estimated using the value of a. Numerically, we find 
AFcRF = (-10±10%)a = -1.4±1.5Wm-^ for the Eocene 
and AFcRF = (70±40%)a = 10.0 ± 6.0 Wm^^ ^j. ^j^g q^.^ 
taceous (if we limit ourselves to p between 1/2 and 1). The 
new estimates are A = O.Se^g jj KW~^m^ for the Creta- 
ceous, and A = 0.42l° ;i^°KW"^m^ for the Eocene. (More 
detail at fig. 3 and table 1). 

3.5. Warming since the last glacial maximum 

Several studies have attempted to estimate the global sen- 
sitivity by comparing the temperature increase since the 
last glacial maximum (LGM) with the radiative forcing 
change. For example, Hoffert and Covey [1992] estimate 
that a radiative forcing of AF — 6.7 ± 0.9Wm^^ is re- 
ponsible for a temperature increase of AT = 3 ± 0.6°K. 
On the other hand, Hansen et al. [1993] find a higher sen- 
sitivity. This is because they find that a similar radia- 
tive forcing of AT — 7.1 ± l.SWm^'^ is responsible for 
a much larger AT = 5 ± 1°K increase. The main differ- 
ence is that Hoffert & Covey base their temperature csi- 
mate on the oceanic CLIMAP temperature reconstruction, 
while Hansen et al. based theirs on land temperature prox- 
ies. We will adopt the average temperature change and 
increase the error to be conservative. Namely, we choose 
AT = 4 ± 1.5°K. Similarly we take AF = 6.9 ± l.SWm'^. 
This gives A = CSSl^ j^^KW-^m^ (as detailed in fig. 2 and 
table 1). 

Again, the above estimates do not include the net ra- 
diative forcing change due to CRF modulation of the cloud 
cover. 

Christl et al. [1998] and Frank et al. [1997] assumed that 
^''Be flux modulation on this time scale is primarily a result 
of modulation by the varying geomagnetic field. Using this 
flux they derived that the geomagnetic field was about 50% 
its present value at 20 ka before present. Moreover, the to- 
tal effect of the terrestrial field is a ~ 20% reduction of the 
CRF penetrating down to ~ 3km, where low altitude clouds 
typically form [Compton, 1933]. Namely, the 50% reduction 
in the global magnetic field should correspond to a reduc- 
tion of about 10% in the relevant CRF. i.e., AFcrf « 0.1a. 
This should be compared with the 7% variation over the 
solar cycle. 

Sharma [2002] relaxed the assumption that the ^"Be flux 
modulation is predominantly terrestrial. By using indepen- 
dent proxies for the terrestrial field, he obtained that the 
field was only 30% lower than today, corresponding to a 
~ 6% reduction in the high energy CRF. The rest of the 
^''Be flux variations, were attributed to a reduced solar mod- 
ulation factor (f) [Masarik and Beer, 1999], that at 20 ka was 
about a 1/3 of its current day average value of ~550 MeV. 
If we take the expression for the CR differential number 
flux J[Kg,<f>) given in Masarik and Beer [1999] and inte- 
grate EpJ{Ep, 4>)dEp to get the total flux reaching the 
troposphere (i.e., with Ep^Eo w IQGeV), we find that the 
increased solar activity as borne out from the increased so- 
lar modulation factor was responsible for an additional 6% 



reduction in the high energy CRF. We therefore get a total 
reduction of about 12% in the CRF, and AFcrf « 0.12a. 
We take this value. 

Since we flnd that a larger total radiative forcing is re- 
sponsible for the same temperature change, we obtain a 
lower estimate, A = 0.48lQ;ig°KW~^ m^ (also detailed in 
fig. 3 and table 1). 

Instead of using the radiative forcing through cloud cover 
modification, we can use our limits of /u which are not based 
on LACC forcing, but instead on the observed temperature 
change over the phanerozoic. We found /tt = 6.5 ± 2.5°K. 
Thus, the 12% decrease in the CRF causing low altitude 
tropospheric ionization should translate into a warming con- 
tribution of ATcRF ~ 4.8 to 10.8°K warming. We now 
use an equation of the form AT = A(ATo -I- AFcrf) = 
AAFo + ATcRF, and obtain: A = 0.47lo:i9°K W"^ m^ . 

To summarize, the effect of the decreased CRF since the 
LGM is to reduce our estimate for A by about 20%, which 
is smaller than the error in the estimate itself. This result is 
valid also if we don't believe the CRF-climate link is through 
LACC moduation, but merely that such a link exists. 

3.6. Warming over the past century 

Climate sensitivity can also be estimated using the global 
warming observed over the past century once the radiative 
forcing with their uncertainties are estimated. 

Since the time scale is relatively short, it is necessary to 
consider the finite heat capacity of the oceans. We base 
our analysis here on the work of Gregory et al. [2002] who 
tackled this problem by considering the heat flux into the 
ocean in the energy budget. The main difference between 
our modified analysis here and that of Gregory et al. [2002], 
is that we will also consider the radiative forcing associated 
with the decreased CRF over the past century. Unlike the 
warming since the LGM, where this was a small correction, 
here it is will prove to be a notable one. 
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Figure 3. The Probability Distribution function for A (and 
ATx2) obtained by comparing radiation budget differences 
to temperature cinange over various time scales, assuming 
that CRF variations do affect the global climate through 
modulation of the low altitude cloud cover, and that the 
relation between cloud cover changes and radiative forcing 
are given by the nominal range. The cases are the same as 
in figure 2 with the added CRF effect. In addition, we have 
case (7) Global sensitivity from correlation between CRF 
variations and temperature variations over the Phanerozoic. 
No equivalent for this case exists in the previous set because 
the CRF/climate link is required to explain this data. And 
(8) Warming since the LGM assuming n is given by the 
CRF/temperature correlation over the Phanerozoic and not 
necessarily through LACC variations. 
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Again, we assume that the CRF modulates the LACC 
and that its radiative forcing is given in §2. 

Gregory et al. [2002] compared the period 1850-1900 with 
1950-1990. Since the CRF record does not go back far 
enough, we need to use proxy data. A good choice is the 
geomagnetic aa activity index. The advantage of the aa in- 
dex is that it can directly compare the long secular trends 
to variations over the solar cycle. In particular, we find that 
the secular increase between the two periods is roughly the 
same as the total variation in the aa index over the solar 
cycle. In both cases it is roughly 12nT. If we further assume 
that the ratio between the secular CRF variations to those 
over the solar cycle is similar to the ratio in the aa index, 
wo find that the decrease in the ~ 10 GcV CRF between 
the above periods should bo about 7%, the typical min to 
max variations over the solar cycle. Using our estimate for 
the radiative forcing change in §2, this corresponds to about 
1.0±0.4Wm-^ 

According to Gregory et al. [2002] , 



AAT = F-Q, 



(18) 



where F is the change in the radiative forcing (i.e., in the 
energy balance) , while Q is the average not heat flux which 
entered the oceans between the two_periods. 

In our modified case F = Fo + Fcrf , where Fo is the 
"standard" radiativo_ forcing that was estimated by Gregory 
et al. [2002] to bo Fo = -0.3 to IWm"^. It includes an- 
thropogenic, volcanic, solar luminosity and aerosol contri- 
butions (with the last one contributing the largest uncer- 
tainty). Q was estimated to be 0.32 ± 0.15Wm~^ while 
AT = 0.335 ± 0.033° K (all at 2ct). 

Like Gregory et al. [2002] , we assume that the errors have 
a Gaussian distribution. Following their procedure, we cal- 
culate the probability distribution function (PDF) for A. 
They are given in figs. 2,3, for the CRF and no-CRF cases. 
The added complication i_n the CRF case is the extra PDF 
for a, which implies that Fcrf has a PDF itself. The PDF 
obtained for the Fcrf = case is the same as the result of 
Gregory et al. [2002]. 

Inspection of fig. 3 and table 1 reveals that A = 
0.32lo;io°KW~^m^ at 1-cr confidence. This is a clear re- 
duction from the results of Gregory et al. [2002], where the 
lower 16th percentile for A is 0.67°K m^ and there is no 
formal upper limit (This assumes a prior that A cannot be 
negative) . 

3.7. Variations over the solar cycle 

The shortest time scale we study is that of the solar cy- 
cle. Since the expected signal arising from solar variability 
is small (~ 0.1°K) looking at the recent few cycles is prob- 
lematic since internal variations (such as volcanic eruptions, 
ENSO and other oscillations as well as simple inter-annual 
variations, or "internal noise") in the climate are large and 
can drown the solar signal. To overcome this problem, we 
will look at a much longer temperature record. In particular, 
we use the post-Little Ice ago data (i.e., last 300 yrs) of the 
Jones et al. [1998] thousand year long reconstruction of tem- 
perature for both hemispheres, which includes proxy data of 
tree rings, ice cores, corals, and historical documents. The 
catch is that the solar cycle is not stable, and the actual pe- 
riod varies between about 9 and 12 years. To overcome this 
problem, we do not perform a harmonic analysis. Instead, 
we fold and average the data over the varying solar cycle 
period. Namely, each year is asigned a phase 4> defined as 
the time since last sunspot minimum devided by the length 
of the particular cycle encompassing the given year. Once 
we do so, we can average all points within a phase bin and 
obtain the average temperature. Also, the internal variance 
in the temperature allows us to estimate the error in AT(0). 
The result is depicted in figure 4. 



Evident from the figure is the fact that the temperature 
has a near sinusoidal behavior. By performing a fit to a 
form AT(<;/>) = AT -|- (a/2) cos(27r((^-0o)) with 00 the phase 
relative to the occurrence of maximum sunspot number, we 
find at the 1 — a confidence limit that: 



a = 0.09±0.03°K, 
po = 1.0±0.6rad. 



(19) 
(20) 



The value of 0o implies that the average temperature lags 
behind the maximum sunspot number by 1.8 ± 1.0 years. 
This is to be expected, because of Earth's finite heat capac- 
ity. As a result, the response to radiative perturbations is 
not only damped, but it is lagged as well. 

Other analyses estimated the surface temperature varia- 
tion over the 11-yr solar cycle. Douglass and Clader [2002] 
found a = 0.11 ± 0.02°K, while Whtte et al. [1997] found 
a = 0.10 ± 0.02°K. Together with the current result, we 
will adopt a = 0.10 ± 0.02° K for the temperature variation 
between solar minimum and solar maximum. 

We now use the results of §2. In particular, the above 
temperature variations are assumed to arise from the ~ 1.6% 
variations observed in the LACC, such that the forcing over 
the solar cycle is AFlacc = 1.0±0.4Wm~^. An additional 
contribution of AFfiux ~ 0.35Wm~^ is due to changes in 
the solar fiux. The sensitivity itself is then given by: 



A = 



AT/d 
AF 



(21) 



where d is a damping factor which arises from the finite heat 
capacity of the cfimate system and its inability to reach equi- 
librium at a finite time. 

The value of the damping factor is not well known. In 
principle, it can be obtained in climate models, but these 
give a range of values. Using a simple Ocean/Climate 
model, Schlesinger et al. [2004] obtained a damping of about 
0.25 on the 11-yr solar cycle time scale (and about 0.75 on 
the centennial time scale). Other more extensive simula- 
tions find that the 11-yr solar cycle is damped relative to 
the centennial scale by a factor of ~ 0.54 [Cubasch et al., 
1997], 0.33 [Rind et al., 1999] or by comparing solar forc- 
ing to actual climate responce, to ~ 0.68 [Waple et al., 
2002] . If we further consider that the centennial time scale 
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Figure 4. Plotted are average number of sunspots, the 
global temperature and a sinosoidal fit. The global surface 
temperature is the reconstructed temperature over the past 
300 yrs (with variations on time scales longer than 30 yrs 
removed) folded over the solar cycle length. The average 
sunspot number folded over the same cycle. 
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Table 1. Limits on the sensitivity of global temperature to radiative forcing, using different methods while 
assuming that CRF does or does not affect climate. ATx2 is calulcated using a forcing of 3.71 Wm^^ [Myhre 
et al., 1998]. The columns denote the values below which A (or ATx2) are expected to be found at the given 
probability, as obtained with the particular method, or, once the methods are combined. Either by assuming 
A is strictly constant, or by assuming it ca 

Period (Method) A[°KW-im^] ATx2[°K] 

1% 16% 50% 84% 99% 1% 16% 50% 84% 99% 



Without the effect of Cosmic Rays 



Phanerozoic (CO2/T) 


0.02 


0.19 


0.36 


0.63 


1.05 


0.1 


0.5 


1.3 


2.3 


3.9 


Cretaceous" 


0.29 


0.38 


0.57 


1.01 




1.1 


1.4 


2.1 


3.8 


* 


Eocene'' 


0.03 


0.21 


0.37 


0.56 


0.87 


0.1 


0.8 


1.4 


2.1 


3.3 


LGM''''= 


0.10 


0.38 


0.58 


0.87 


1.48 


0.4 


1.1 


2.2 


3.2 


5.5 


20* Century'* 


0.34 


0.67 


1.31 


* 


* 


1.4 


2.5 


4.9 


* 


* 


iDUltLl v^ycic 


24 


0.58 


0.94 


1.85 




n Q 


2 2 


0.0 


u.o 


* 


Conibiuod (A = const) 


0.33 


0.44 


0.52 


0.62 


0.79 


1.3 


1.6 


1.9 


2.3 


2.9 


Combined (X = Ao + hAT) 


0.21 


0. 1:-! 


0.51 


().()() 


0.S7 


0.9 




2.0 


2.r, 


:-!.2 




With 


the effect of 


Cosmic Rays 












Phanerozoic (C02/T)'= 


0.02 


0.14 


0.28 


0.43 


0.65 


0.1 


0.5 


1.0 


1.6 


2.4 


Phanerozoic (CRF/T)-'' 


0.02 


0.22 


0.44 


0.89 


* 


0.1 


0.8 


1.6 


3.3 


* 


Phanerozoic (Clouds)^'-' 


0.22 










0.8 










Cretaceous"'-'^ 


0.13 


0.24 


0.36 


0.60 


1.94 


0.5 


0.9 


1.3 


2.2 


7.2 


Eocene'''-''''' 


0.03 


0.24 


0.42 


0.67 


1.32 


0.1 


0.9 


1.6 


2.5 


4.9 


LCM (LACC)''''^-'' 


0.08 


0.30 


0.47 


0.69 


1.12 


0.3 


1.1 


1.7 


2.5 


4.1 


LGM (//)'''='9 


0.04 


0.28 


0.47 


0.74 


1.30 


0.2 


1.0 


1.8 


2.8 


4.8 


20*'' Century-'' 


0.16 


0.22 


0.32 


0.64 


* 


0.6 


0.8 


1.2 


2.4 


* 


Solar Cycles 


0.08 


0.15 


0.26 


0.52 


2.30 


0.3 


0.6 


1.0 


1.9 


8.6 


Combined'^A = const) 


0.17 


0.24 


0.29 


0.37 


0.50 


0.6 


0.9 


1.1 


1.4 


1.9 


Combined''(A = Ao + bAT) 


0.25 


0.27 


0.35 


0.45 


0.68 


0.9 


1.0 


1.3 


1.7 


2.5 



" Based on Hoffert and Covey [1992]. '' Based on Covey et al. [1996]. " Based on Hansen et al. [1993]. 

Based on Gregory et al. [2002]. " Based on Shaviv and Veizer [2003]. 
f Assumes tlic CRF-climate is tlirough modulaton of the LACC, witli a = 14 ± 6Wm~^ 
s Using ATcRF and ^ from the Phanerozoic data instead of AFqrf and a. 

^ The combined PDF does not include the LGM (/i) estimate or the Phanerozoic/Clouds bounds. 
' Limits for ATx2 larger than 10° K are meaningless and therefore not quoted. 

^ The lower limit obtained from the maximum cloud cover changes depends on systematic errors. The confidence limits are 
therefore meaningless. 



is damped relative to the long term response, by a factor of 
~ 0.7-0.75 [Schlesinger et al, 2004; IPCC, 2001], then the 
4 estimates for the damping factor are encompassed within 
d = 0.35±0.15, for periodic oscillations with an 11-yr period. 

Note that by resorting to GCMs for the estimate of the 
damping factor, we are somewhat unfaithful to the spirit of 
this work, which is to estimate the sensitivity independently 
to the usage of GCM simulations. Nevertheless, the damp- 
ing we use is a characteristic describing the relative behavior 
of different time scales. We still avoid using the absolute sen- 
sitivity obtained in GCMs. Moreover, the analysis of [ Waple 
et al, 2002] does indicate that empirically, the damping fac- 
tor is consistent with that obtained in GCMs. The fact that 
this result is somewhat larger than GCMs on average, would 
imply that we maybe underestimating the damping factor, 
and with it, overestimating the climate sensitivity. 

For nominal values, wo find that A = 0.26to ^^"KW"^ m^. 
Without the effects of the CRF, a much larger sensitivity 
(of A = 0.94l",;3p,°KW"-^m^ ) is obtained because the same 
temperature variations arc then to be explained only by the 
relatively small solar flux variations. 

3.8. Combined Results 

We now proceed to combine the PDFs obtained in the 
cases described in fig. 2, when the CRF/climate link is 
neglected, and the cases described in fig. 3, when the 
CRF/LACC effect is included. We combine the results in 
two cases. In the first, we assume that the global tempera- 
ture sensitivity is constant, namely, that it does not depend 
on the average terrestrial temperature. In the second case, 
we allow the sensitivity to be temperature dependent. 

3.8.1. Constant Sensitivity 



When combining the results under the assumption that 
CRF does not introduce a radiative forcing, we can simply 
multiply the PDFs and renormalize the result. The reason 
is that the error in the estimates of all the AFq and AT are 
presumably uncorrclated with each other, and also because 
wo have no prior on the value of A (except perhaps that it 
should be positive). 

On the other hand, when combining the cases which in- 
clude the CRF/LACC effect, we must bear in mind that 
some of the error arises from the uncertainty in a - the rela- 
tion between cloud cover changes and radiative forcing. This 
uncertainty enters 5 PDFs, and we cannot simply multiply 
them. To overcome this obstacle, we calculate the PDFs 
assuming a given a. Then, the combined PDF is given by 



PalliX) ^ 



/[nti-P'(^.Q)]-Pa(«)^Q 



(22) 



Again, this assumes that we have no prior on A, and that 
besides the dependence on a, the PDFs are not in anyway 
correlated with each other. 

In figs. 6,7 and table 1, we plot and describe the com- 
bined PDFs obtained in the two cases. We find that A = 



0.52 



+0.10° 

-0.08 



KW^^m^ if the CRF/Climate fink is neglected, 
and that A = 0.29l";"g°KW"^ m^ if the CRF/LACC link is 
included. Values of upper and lower limits on A and ATx2 
at different confidence limits are given in Table 1. 
3.8.2. Variable Sensitivity 

We now alleviate the assumption of a constant sensitivity 
and allow a linear relation in the form A = Ao -|- bAT, where 
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Figure 5. The estimated sensitivity A as a function of the average temperature AT relative to today over 
which the sensitivity was calculated (e.g., average temperature between today and a given epoch if conditions 
at the given epoch and today are used to estimate the sensitivity). The values are for the Last Glacial Maxi- 
mum (LGM), 11 year solar cycle over the past 200 years (11), 20*'' century global warming (20), Phanerozoic 
though comparison of the tropical temperature to CRF variations (Phi) or to CO2 variations (Ph2), Eocene 
(Eo) and Mid-Cretaceous (Cr). Panel (a) assumes that the CRF contributes no radiative forcing while panel 
(b) assumes that the CRF does affect climate. Thus, the "Phi" measurement is not applicable and does not 
appear in panel (a). From the figures it is evident that: (i) The expectation value for A is lower if CRF affects 
climate, (ii) The values obtained using different paleoclimatic data are notably more consistent with each 
other if CRF does affect climate (iii) There is no significant trend in A vs. AT. 



Ao is the sensitivity today and b = d\/dT. Because the er- 
rors are not Gaussian (the distributions are generally skewed 
towards higher A's) we cannot fit A(A) using a simple linear 
least squares. Instead, we calculate 



P u(Xo b) = I[I]^Umo + bAT.,^)]P4a)da 

J [[f^^^P,i\o + bATi,a)]Pc{a)dadXdb 

and Paa(Ao) = / PaiiiXo,b)db. 

Here we find that A = 0.54l°;f "KW-^m^ if 
the CRF/Climate link is neglected, and that A = 
0.35i;];J'^°KW-^ if the CRF/LACC link is included. Val- 
ues of upper and lower limits on A and ATx2 at different 
confidence limits are given in Table 1. 

This is our best estimate for the global sensitivity. It 
translates into a CO2 doubling temperature change of 
Arx2 = 1.3±0.4°K. With the CRF/climate effect neglected, 
this rmmber is ATx2 = 2.0 ± 0.5°K. 

Another point worth mentioning is the fact that once the 
CRF/LACC climate link is included, the median values for 
A obtained using different periods differ from each other by 
typically 1-a or less, while without the CRF/climate con- 
sidered, differences can be larger than 2-a. Namely, the 
CRF/climate effect markably improves the consistency of 
the data. This can be seen in fig. 5. 

4. Discussion 

We compared the radiative forcing and temperature 
change over several different time scales, while taking into 
consideration the alleged link between CRF variations and 
temperature change. We found that the 6 different time 
scales can be used to place similar bound on the global cli- 
mate sensitivity and, when possible, also on the quantitative 
relation between CRF variations and temperature changed. 

Different time scales and methods suffer from different 
uncertainties, which have to be kept in mind. These include: 



1. Limits based on the geocliemical record over the past 
550 million years implicitly assume that estimates of temper- 
ature variations using 6^^0 is well known. In principle, var- 
ious biases might distort this relation and produce a wrong 
temperature scale. 

2. Limits based on large temperature variations and ra- 
diative forcing, assume that one can deduce the sensitiv- 
ity today, to small perturbations, from temperature changes 
associated with large radiative forcings, or worse, when dif- 
ferent conditions existed under which the sensitivity could 
have been different (e.g., over the past 550 milUion years 
ice sheets where mostly absent). To address this point, we 
allowed for A to be a linear function of temperature, but 
found no significant dependence. 

3. Most of the estimates of the climate sensitivity assume 
the CRF/climate link is through modification of the LACC 
and that the radiative forcing associated with it is known. 
This entails in it several assumptions: (a) That LACC 
changes observed by the ISCCP measurements, for exam- 
ple, are indeed well represented by changes in the amount 
of cloud cover as opposed to other cloud characteristics. (2) 
That the incremental cloud cover changes behave as the av- 
erage. (3) That tropospheric ionization does not markably 
affect the global temperature through an effect not related 
to cloud cover modifications. These uncertainties will not 
be resolved without detailed understanding of how and to 
what extent does atmospheric ionization affects the format 
tion of cloud condensation nuclei, which affect cloud cover. 
Nevertheless, the value of the sensitivity could still be in- 
dependently bracketed using argumentation which does not 
assume the relation is through the LACC variations. 

4. Over the short time scales associated with the solax 
cycle, one of the main uncertainties is the damping factor in 
the effect of a changed radiative forcing. 

Having said that, the fact that about half a dozen in- 
dependent analyses based on paleoclimatic to recent data 
yield roughly the same sensitivity should indicate that we 
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are probably not missing large radiative forcing terms. Oth- 
erwise, there is no reason, other than chance, to obtain re- 
sults consistent with each other. Moreover, the notably im- 
proved agreement between the sensitivities obtained once 
the CRF/climate effect is taken into account, is yet another 
indicator that the effect is real (see panel A vs. B in fig. 5). 

Our best estimate is ATx2 ~ 1.3°K (0.9°K < ATx2 < 
2.5°K at the 99% confidence level). This is at the lower end 
of the often quoted range of ATx2 = 1.5 to 4.5°K [IPCC, 
2001] obtained from Global Circulation Models (GCMs). 
Cess et al. [1989] have shown that the climate sensitivity 
obtained in this type of simulations predominantly depends 
on how clouds are treated, and whether they contribute a 
positive or negative feedback. The models roughly give that 
« 2.2Wm-^ /°K- AF^ioud/Af with AFdoud being the 
feedback forcing of clouds associated with a temperature 
change of AT. Thus, for a GCM to be compatible with the 
results obtained here, a negative cloud feedback is required. 
One such example was suggested by Lindzen et al. [2001] 

On the other side of the coin, we can also rule out very 
small climate sensitivities. This can be used for example 
to place a limit on possible large negative feedbacks, or to 
a lower limit on the effect of anthropogenic greenhouse gas 
(GHG) warming. 

Since the beginning of the industrial era (~ 1750), non- 
solar sources contributed a net forcing of 0.85 ± 1.3Wm~^ 
[IPCC, 2001] (assuming the errors are Gaussian). Over the 
past century alone, this number is 0.5 ± 1.3Wm~^. The 
main reason why the error is large is because of the uncer- 
tain "indirect" contribution of aerosols, namely, their effect 
on cloud cover. It is currently estimated to be in the range 
-1 ± IWm"^ [IPCC, 2001]. Thus, anthropogenic sources 
alone contributed to a warming of 0.13 ± 0.33°K since the 
beginning of the 20**^ Century. 

The sensitivity result can also be used to estimate the 
solar contribution towards global warming. Over the past 
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Figure 6. The combined probability distribution functions 
for A obtained by combining the PDFs given in fig. 2 and 
fig. 3 when the CRF/climate link is either neglected or in- 
cluded. (In the latter case, the combination is done as ex- 
plained in the text). Thin lines denote the result if A is 
assumed to be temperature insensitive, while the heavy line 
is the result obtained when A is allowed to be a linear func- 
tion of the temperature. The the latter case, the distribution 
for A today is plotted. Also marked are the two additional 
constraints obtained from the Phanerozoic data which do 
not depend on the CRF/LACC link, as well as the sensi- 
tivity range of 1.5 to 4.5°K, which according to the IPCC 
[2001] is "widely cited". Note that the sensitivity range of 
the 15 GCM models actually used by the IPCC [2001] is 2.0 
to 5.1°K. 



century, the increased solar activity has been responsible for 
a stronger solar wind and a lower CRF. Using results of §3.6, 
the reduced ionization and LACC were responsible for an 
increased radiative forcing of 1.0 ± 0.4Wm~^. In addition, 
the globally averaged solar luminosity increased by about 
0.4±0.1Wm~^ according to Solanki and Fligge [1998]; Hoyt 
and Schatten [1993]; Lean et al. [1995]. Thus, increased solar 
activity is responsible for a total increase of 1.4±0.4 Wm~^ . 
Using our estimate for A, we find ATsoiar ~ 0.37 ± 0.13°K. 

We therefore find that the combined solar and anthro- 
pogenic sources were responsible for an increase of 0.50 ± 
0.35°K. This should be compared with the observed 0.57 ± 
0.17°K increase in global surface temperature [IPCC, 2001]. 
In other words, the result we find for the sensitivity and 
drivers are consistent with the observed temperature in- 
crease. The solar and anthropogenic contributions are com- 
parable, nevertheless, it is more likely that the former con- 
tribution is somewhat larger, contributing about 2/3s of the 
observed temperature increase. This conclusion was inde- 
pendently reached by comparing the non-monotonic tem- 
perature increase with the non-monotonic solar activity in- 
crease and the monotonic increase in GHGs [Soon et al., 
1996b]. 

We therefore conclude that both GHGs and solar vari- 
ability are the main drivers of global warming, through a 
temperature sensitivity that is notably lower than the whole 
range obtained in GCMs. Its value is surprisingly close to 
the black body value of 0.30°KW'^m^ 
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